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Abstract 

We investigate the segregation of the extragalactic population via colour criteria 
to produce an efficient and inexpensive methodology to select specific source popu- 
lations as a function of far-infrared flux. Combining galaxy evolution scenarios and a 
detailed spectral library of galaxies, we produce simulated catalogues incorporating 
segregation of the extragalactic population into component types (Normal, star- 
forming, AGN) via color cuts. As a practical application we apply our criteria to 
the deepest survey to be undertaken in the far-infrared with the AKARI (formerly 
ASTRO-F) satellite. Using the far-infrared wavebands of the Far-Infrared Surveyor 
(FIS, one of the focal-plane instruments on AKARI) we successfully segregate the 
normal, starburst and ULIRG populations. We also show that with additional MIR 
imaging from AKARI's Infrared Camera (IRC), significant contamination and/or 
degeneracy can be further decreased and show a particular example of the separa- 
tion of cool normal galaxies and cold ULIRG sources. We conclude that our criteria 
provide an efficient means of selecting source populations (including rare luminous 
objects) and produce colour-segregated source counts without the requirement of 
time intensive ground-based follow up to differentiate between the general galaxy 
population. 

Key words: Space-based infrared telescopes, Far infrared, Galaxy evolution, 
Extragalactic surveys 
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1 Introduction 



Fueled by rapid advances in detector technology, enabling fast and accu- 
rate mapping of large areas of the sky, extragalactic astrophysics is enjoy- 
ing a golden age where astronomers have become spoiled for choice with 
the enormous statistical survey data sets available probing out to cosmo- 
logically significant distances such as the optical Sloan Digital Sky Surve y 



(SDSS) covering one quarter of the sky (lAdelman-McCarthy et al.l . 120061 ) 



the near infrared UKIRT Infra red Deep Sky Survey (UKIDSS, 4000 square 
degrees ([Lawrence et al. 1. 120061)1 and the Spitz er Wide-area Infrared Extra- 
galactic Survey (SWIRE. [Lonsdale et al.l tooi )). The most basic tool for the 
study of cosmological galaxy evolution is the monochromatic (single band) 
number count of extragalactic sources. However aside from source counts and 
clustering analysis of large scale structure, deeper investigation into the extra- 
galactic population and cosmological evolution will inevitable require multi- 
wavelength analysis and hence galaxy colours. Ultimately we would like to 
define the extragalactic zoo species by species and derive their corresponding 
distances to determine the star-formation rate and sequence of structure for- 
mation in the Universe. However such studies often involve extremely time 
intensive follow up campaigns with either deep imaging or spectroscopy. As 
an intermediate step, galaxy colours have been frequently used to gain a 
cheap yet effective insight into galaxy populations and rough distances at 
a fraction of the cost of large scale follow-up campaigns. Such colour anal- 
ysis has been extremely effective in selecting and segregating populations of 
sources in the near-infr ared, optical and u l tra-violet waveband s (e.g. the Ly- 



man bre ak population ( Steidel et al. . 1996 ; Madau et al. . 19961 ). ERO & BzK 



sources ( lElston et al.l . Il988l ; ICimatti et al.l . 120021 )). At infrared wavelengths 
astronomers have often had to work with fewer bands and less well behaved 
detectors, yet the colour criteria has also been applied with great success 
i n the mid-far-infr ared with the pioneering I RAS mission's all-sky survey 
( Soifer et al.l.ll987l). U sing simple colour criteria. iRowan- Robinson fc Crawford 
( 119891 ). iHelou I (119861 ) created a parameterization of the extragalactic popu- 
lation based on the IRAS infrared colours in four bands. Within this frame- 
work, a normal galaxy component was defined by cool lOO/mi/60/im colours, 
a starburst component defined by warm lOO/mi/60/mi colours extending to 
an ultraluminous component at higher infrared luminosities and an AGN 
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component was denned by high 25 / um/60/im colours due to strong emission 
from the mid-infrared dust torus. Two decades after the IRAS all-sky sur- 
vey, the Japanese AKARI mission (launched successfully on Februrary 21st 
2006, formerly known as ASTRO-F) will perform a new All-Sky Survey in 
4 far-infrared and 2 mid-infrared bands to much improved sensitivity, higher 



spatial resolution and wider wavel ength coverage (IShibail . 12004 iPearson et al. 



2004 ; iPearson. Shibai et all 120070 . AKARI is equipped with a 68.5 cm cooled 
telescope (IKaneda et all 120051 ) and two focal plane instruments, the Far- 
Infrared Surveyor (FIS) and the InfraRed Camera (IRC). The FIS has two 
2-dimensional detector arrays and observes in four far-infrared bands between 
50 and 180 /im ( centred on N6 (65u m), WIDE-S (90/iin), WIDE-L (140/rni), 
N160 Q60/W). Kawada et al.l (120040 ) . The IRC consists of three cameras 
covering 1.7- 26 fim in 9 bands (N2 (2.4/iin), N3 (3.2/um), N4 (4.1//m), S7 
(7/um), S9W (9/iin), Sll (11/mi), L15 (15/im) L18W (18/W), L2 4 (24/im)) 
with fields of view of approximately 10' x 10' ( Onaka et all 2004 ). Both the 
FIS and IRC instruments also have capabilities for spectroscopy. In addition 
to the All-Sky Survey, AKARI will also perform deep pointed surveys over 
areas of one t o ten square degrees at the North (NEP) and South (SEP) 



Ecliptic Poles (IMatsuhara et al.l . 120061 ). These major legacy surveys are ex- 



pected to yield hundreds of thousands of sources and thus require a major 
data reduction effort. As a precursor to these surveys, the aim of this work is 
to investigate the feasibility of colour segregation of the infrared populations 
with the minimum peripheral data and effort as possible in a similar manner 
to the IRAS survey. Our colour segregation criteria will have merits not only 
for the AKARI data sets but also for similar large data sets expected from the 
Spitzer & Herschel missions. In this paper we concentrate on the deep region 
selected for the AKARI SEP survey while a more d etailed investigati o n of o ur 
methods and their applications will be reported in lPearson &: Jeongi (120071 ). 



2 Phenomenological Model for Colour Segregations 



2. 1 Spectral Energy Distribution Library 



To model the colours of the extragalactic populations we have gathered a selec- 
tion of spectral energy distribution (SED) templates for various galaxy popula- 
tions from four large readily available contemporary spectral libraries. Compu- 
tational constraints from our simulations restricts us from utilizing entire sets 
of thousands of SEDs so rather we have selected a moderate sub-sample which 
spans the entire range of spectral properties of the libraries. The selected sets 
of SEDs specifically model the quiescent normal galaxy population, the star- 
burst population (increasing in luminosity: starburst, luminous infrared galax- 
ies (LIRG), ultra-luminous infrared galaxies (ULIRG)) and the AGN popula- 
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tion. To model the quiescent normal galaxy po pulation we have selected spec 



Efstathiou fc Rowan-Robinson 



Efstathiou fc Rowan-Robinson 



tral te mpla tes from th e libra ries described in 
(|2003h and bale et all fl200lh . The models of 
( 120031 ) provide good fits to the IR AS, ISO and Spitzer selected sample s of the 

cool, quiescent galaxy population (IRowan-Robinson et all 12004 12005|) as well 

as th e cold component of the sub-millimetre source population (jEfstathiou &: Rowan- Robinson . 



20031 ) an d the galaxy source c ounts from sub-millimetre to near-infrared wave- 



lengths (jPearsonl . l200ll . 120071 ). The model templates of Dale et al. have been 
shown to provide good fits to loc al normal galaxies from both ISO and more re- 
cently Spitzer se l ected samples ( ID ale et all l200ll . 120051 ) , high-z Spitzer sources 



( Appleton et all 120041 ) and the ISO & Spitzer mid-infrared source counts 



( jPearsonl . 120051 ). The starburst, LI RG, ULIRG population SEDs are taken 

from t he ra diative transfer model s of lEfstathiou. Rowan- Robinson fc Siebenmorgen 

(l2000h and iTakagi et all (j2003al fbh . The starburst models of Efstathiou have 

provided good fits to a wide variety of individual objects from starburst (e.g. 

M82), LI RG (e.g.NGC6090, IRAS1445-4343 ) to ULIRG sources (e.g. ARP220, 

MK231), ( Efstathiou. Rowan- Robinson &: Siebenmorgen . 20001 Efstathiou fc Siebenmorgen . 



2005) and also the recent results from the Spitzer SWIRE survey (IRowan-Robinson et al 



20051 ). The models of Takagi et al . prov ide good fits to the SEDs of dusty sub- 



millimetre galaxies (ITakagi et all 120041) . The AGN population is modelled on 



the re cent SED library of lSiebenmorgen et al.l (120041 ) ; ISiebenmorgen fc Efstathiou 
( 120051 ). The SED libraries provide models for both Type 1 and Type 2 AGN by 
the consideration of tapered disc geometries. The selected S EDs are shown in 
Figure 1. Further details of the model selection are provided in lPearson fc Jeong 
(l2007h 



2.2 Simulated Catalogues and Galaxy Colours 



To estimate the expected fluxes and corresponding colours in the AKAR I 
bands we use the contemporary galaxy evolution model of iPearsonl (120071 ). 
This model uses a backward evolution process based on the the IRAS all- 
sky PSCz multi-component luminosity function defined at 60/mi comprising 
of cool normal quiescent galaxies, and a warmer component defined by in- 
frared luminosity as Lm < 10 u Lq starburst galaxies, Lm > 10 11 Lq LIRG 
sources, Lm > 10 12 Lq ULIRG sources and AGN. This model fits the ob- 
served source counts from 2-1200/im from the IRAS, IS O, Spitzer m i ssions 
and the SCUBA/JCMT, MAMBO/ I RAM instruments (iBertin et all ' 



Puget et al.l. 1199a iDole et all 120041: iFazio et all 120041; lAussel et al 



Serieant et al. I . 1200(1 iPapovich et all 12004 iBlain et all l2002t iGreve et al. 



1997 



1999 



20041 ). Our source colour information is simulate d by passing a cata logue of 



point sources through the instrument simulator of I Jeong et al.l (120041 ) . An in- 
put catalogue is c onstructed for t he 5 broad galaxy types outlined above within 
the framework of IPearsonl (120071 ). The catalogue is constructed as a function 
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of increasing redshift in incremental redshift bins. For each redshift bin, the 
theoretical number-redshift distribution from the model is used to randomly 
select the appropriate number of galaxies at that redshift and a corresponding 
luminosity by sampling the input luminosity function. Each source from the 
catalogue is assigned a SED from the library depending on the galaxy type 
and the expected flux for each source from the corresponding SED at the de- 
sired wavelength and redshift is calculated and convolved with the appropriate 
filter response function. A one sigma variation in the SED flux (i.e. for a given 
source in the catalogue the assigned flux at a given wavelength is selected 
randomly between ±lcr of the original SED template flux) is also included at 
this stage to allow for a wider range of output fluxes for a limited number 
of input SEDs. Note that although the i nstrument simulator is optimized for 



the AKARI mission (IJeong et all [2004J), it is a simple matter to extend the 



simulator to encompass the characteristics of t he infrared pass-b a nds o f other 
space missions (e.g. Spitzer, Herschel, etc., see iPearson fc Jeongi (120071 )). The 
final output from the simulator is a catalogue of sources with associated SED, 
and filter band convolved fluxes at a set of desired wavelengths. 



3 Predicted Colours of the Extragalactic Population in the AKARI 
SEP survey 



In this work, the simulated source catalogue is used to investigate the deep 
AKARI SEP survey. The AKARI SEP survey is a guaranteed time AKARI 
Mission Program (MP). The current projected area spans ~8 square degrees in 
a low cirrus brightness region centred at R.A.=4 ,l 44 m 00 s , Dec=-53deg 20'00."0 
(J2000) where the cirrus brightness is spectacularly low < B wo >~ 0.2 MJy/sr 



( ISchlegel. Finkbeiner fc Davis I . Il998l ). This survey is estimated to reach the 



source confusion limit in the far- infrared FIS wide bands a nd will be th e deep- 
est far-infrared image of the Universe to date (note that IJeong et al.l (120061 ) 
predict confusion limits of ~10mJy & 60mJy in the FIS 90/im and 140/im 
bands respectively based on the source confusion criterion of 20 beams per 
source). Shallow mid-infrared data will also be available over the entire area, 
along with deep optical i?-band imaging to 25th magnitude. In addition, an 
area of approximately 1 square degree will be covered by IRC mid-infrared 
imaging to ~100/xJy. Source confusion at optical wavelengths may also pose 
complications for the identification of faint infrared sources. In the optical 
R-band, the source density may be up to an order of magnitude higher than 
the far-infrared AKARI source density. Furthermore, although the advent of 
wide field optical multi-object spec troscopy inst r umen ts, such as AAOmega 
on the Anglo- Australian Telescope (ISharp et al.l . 120061 ) will make large scale 
follow up a reality, note that at these optical magnitudes (B<23) there will 
still be ~5 optical sources per AKARI beam. Therefore, any identification / 
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population segregation that is possible using the infrared data alone will be 
of great value. 



In Figure 2 the results of our analysis are presented. The top panel shows 
th e predict e d sou rce counts over lOmJy-Uy & 50mJy-lJy from the models 



of iPearsonl (120071 ) in the AKARI FIS 90/on (left) & 140/im (right) bands 
respectively. The AKARI SEP survey should probe to depths of ~10mJy & 
50mJy (3a) in the 90/im & 140/im bands respectively. In the 90/iin band, 
at brighter fluxes (S >100mJy) the source counts are predominantly domi- 
nated by cool normal galaxies due to the dust emission hump peaking in the 
range ~100-200/im. At fainter fluxes (S <100mJy) the warmer evolving star- 
burst and LIRG populations dominate. The 140/im band, samples the cooler 
Rayleigh- Jeans regime of the dust emission hump and is dominated by cool 
normal galaxies and cold ULIRG sources to S ~100mJy. The fainter fluxes 
again see the emergence of the warmer LIRG population. 

The middle panels of Figure 2 show an example colour criterion utilizing the 
AKARI WIDE bands. The 140/im/90/im and 90/xm/140/xm colours are plot- 
ted as a function of 90/im & 140/xm flux respectively on the same flux scales 
as the source counts in the top panel in the figure. Limits in the flux and 
the colour corresponding to a 3a 90/iin flux limit of lOmJy and a 140/im flux 
limit of 50mJy are plotted as dotted lines. Immediately it can be seen that 
the far-infrared colours provide a reasonable discriminator of galaxy type with 
the more active star-forming galaxies occupying the parameter space corre- 
sponding to warmer colours (higher 90/im flux correspond ing to 140/xm/90/ar i 
colours < 1.5 and 90/xm/140/xm colours > 1). Note that Frayer et al.l ( 20061 ) 



have carried out a similar exercise with sources from the Spitzer First Look 
Survey (FLS) although their approach was instead to fit different dust tem- 
peratures to the far-infrared colours using the MIPS 160/mi and 70/tm bands, 
interpreting a cooler dust temperature as cold infrared colours. AKARI, with 
its 4 FIS bands will have a greater power than Spitzer in such far-infrared 
population colour segregation. In addition to providing a crude separation of 
cool normal galaxies and warmer starburst and LIRG sources, the FIS colour 
segregation is also capable of identifying a distinct cold ULIRG population. 
These sources occupy a unique colour-flux parameter space in Figure 2. For 
the case of the 140/im band (right-hand middle panel), the ULIRGs occupy 
the parameter space defined by S(140/xm)>100mJy with 90/mi/140/mi colour 
<0.2. In the case of the 90/xm band the ULIRGs are found at faint er fluxes 



of S (9 0/mi) < 30mJy with 140/mi/90/im colours > 6. The models of [Pearson 



( 120071 ) predict that in the 90/im band, these colour selected ULIRGS will be 
high redshift (l<z<3) sources whilst the 140/im colour selected ULIRGs will 
be predominantly at redshifts z<1.5. The fainter 90/im ULIRG sample should 
include almost all the 140/im ULIRGs. From the middle panel in Figure 2 it 
can also be seen that the ULIRG samples selected in both the 90/im & 140/mi 
bands descend and ascend respectively into the colour parameter space occu- 



6 



pied by the normal galaxies, i.e. that there is contamination of the cold ULIRG 
population by the cooler colours in the normal galaxy population. This con- 
tamination of ULIRGs by cool (low redshift) normal galaxies has also been 
reported in other colour selection studies over far-mid-infrared wavelengths 
( ITakagi fc Pearson! . 120051 ). Using suitable additional colour criteria it is also 
possible to de-convolve the ULIRG sample from the cool normal galaxy pop- 
ulation. The bottom panel of Figure 2 gives an example of the potential power 
of additional colour constraints to further segregate the galaxy populations. 
The bottom panel shows the segregation of the normal and ULIRG popula- 
tions when one of the deep mid-infrared bands is used (in this case the IRC S7 
7/iin band) in conjunction with the FIS bands. The bottom left panel shows 
the results for the 7yum/90/im colours as a function of 90/zm flux. The nor- 
mal galaxies are separated from the cold ULIRG population by virtue of their 
higher mid- infrared/far- infrared flux ratio (see Figure 1). The slight remaining 
contamination that can be seen is actually from a hotter ULIRG population 
whose mid-infrared emission is dominated by an AGN. An even more signifi- 
cant segregation can be achieved using the 7^m/140^m colours as a function 
of the longer wavelength 140/mi flux. Since the far-infrared emission of the 
cold ULIRGs is redshifted into the 140//m band, the mid-infrared/far-infrared 
flux ratio is even lower. In addition, the normal galaxies exhibit a smaller 
scatter in 7 / um/140/im colours as a function of 140/um flux since there is less 
variation in the SEDs on the Rayleigh- Jeans slope than at the shorter wave- 
lengths sampled by the 90/im band. The bottom panels of Figure 2 also show 
that, although to a lesser extent, a similar segregation can also be achieved 
for the AGN population. 



4 Conclusions 



Using a suite of SED libraries combined with a contemporary galaxy evolu- 
tion model and instrument simulator, we have carried out a pilot study of, 
and demonstrated a method to segregate galaxy populations in large data 
sets in a similar manner to the methods utilized for the IRAS all-sky survey. 
Using simple colour criteria, the normal galaxy, starburst/LIRG & ULIRG 
populations can be succesfully discriminated as a function of far-infrared flux 
without the need for large-scale optical follow up campaigns. This method has 
been applied to the AKARI deep survey at the SEP where it is successful 
in achieving a general segregation of the various extragalactic populations. 
The addition of more bands with a greater dispersion in wavelength can be 
used to further refine and decontaminate specific target populations (in this 
example the ULIRG and normal galaxy populations). The colour segregation 
criteria can thus be used to construct basic multi-component source counts 
without the need for expensive ground-based follow up. This is illustrated in 
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Figure 2. For the AKARI FIS 90/iin band, the normal galaxies are identi- 
fied with sources having colours over the range 1.5<140/im/90//m<5 over all 
values of 90/im flux. The starburst/LIRG sources are identified with sources 
of colours 90/im/140/im colours > 1. The ULIRG sources can be identified 
by the dual colour criterea of 140/mi/90/im>6 & 7,um/90yum<0.007. In the 
140/mi band, the normal galaxies can be identified with sources of colour 
0.2<90/xm/140/xm<0.7 while the starburst & LIRG sources have generally 
warmer 90^tm/140^tm colours >0.7. The ULIRGs in the 140/mi band are de- 
fined as sources of colour 90^tm/140^tm<0.2 & 7/mi/140 / um<0.02 over a flux 
range S(140 / um)> lOOmJy. Although not the focus of this present work, a 
similar criterion can also be applied to the AGN population. 

The colour-segregation method and corresponding colour segregated source 
counts holds great promise for the data sets that are expected from the large 
area surveys that are being conducted or will be conducted by the present and 
next generation of space-borne far-infrared telescopes such as Spitzer, AKARI 
& Herschel. A more detailed description of this work and application to the 
A KARI All-Sky Survey an d possible surveys with Herschel will be presented 
in 



Pearson fc Jeongi (120071 ) 
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Fig. 1. Selected spectral energy distributions (SED) from the spectral libraries 
described in the text normalized to an arbitrary flux at 60/iin. (a) Selected normal 
galaxies from the libraries of Efstathiou and Dale (b) Selected starburst galaxies 
from the libraries of Efstathiou and Takagi. (c) Selected LIRG sources from the 
libraries of Efstathiou and Takagi. (d) Selected ULIRG sources from the libraries 
of Efstathiou and Takagi. (e) Selected AGN from the libraries of Siebenmorgen and 
Efstathiou. 
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Fig. 2. Example of colour segregation criteria. Top panel show s the pr e dicted 
source counts over lOmJy-Uy & 50mJy-lJy from the models of iPearsonl (|2007l ) 
in the AKARI FIS WIDE-S 90/jm (left) & WIDE-L 140/im {right) bands respec- 
tively. Results are shown for total counts and component counts (Normal, starburst, 
LIRG, ULIRG, AGN populations). Middle panels show an example colour criterion 
utilizing the AKARI WIDE bands. The 140/im/90/im and 90/zm/140/xm colours 
are plotted as a function of 90/im Sz 140/im flux respectively on the same flux scale 
as the source counts above. Bottom panel shows the segregation of the normal and 
ULIRG/ AGN populations when one of the deep mid-infrared bands is used (IRC 
S7 7fim band) in conjunction with the FIS bands. In the figure, the green lines and 
green open circles represent the normal galaxies, the blue lines and blue crosses the 
starburst galaxies, pink lines and pink plus symbols the LIRG sources, red lines and 
red open circles the ULIRG sources and the black lines and black filled circles the 
AGN. The diagonal black dashed lines in the bottom 4 panels correspond to the 

flux limits for the SEP survey in the relevant AKARI band. 
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